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SUMMARY 


The  research  performed  under  this  contract  was  a  feasibility  study  in  the 
development  of  applications  of  new  microcomputer  graphics  technology  to  the 
analysis  and  interpretation  of  clinical  pharmacological  data.  This  involved 
continuing  development  of  comprehensive  programs  for  analysis*  Interpretation, 
and  simulation  of  pharmacoki neti c  data*  dose-response  kinetic  data*  and  other 
data  relevant  to  new  drug  development*  for  use  with  the  Tektronix  4052 
Microcomputer  Graphics  System.  The  combination  of  such  modern  analytical  and 
illustrative  methods  in  clinical  pharmacology,  based  on  new  high-speed 
microcomputers  and  associated  graphics,  are  thought  to  greatly  reduce  both 
cost  and  time  involved  In  the  overall  process  of  clinical  evaluation  of  new 
drugs  In  the  U.S.  Army  Drug  Development  Program. 

The  work  performed  under  this  contract  includes  the  following: 

1.  Development  of  a  program  for  the  analysis  of  linear  pharmacokinetic  data 
based  on  the  method  of  separate  exponentials  and  using  compulsory  graphics. 

2.  Development  of  a  program  for  the  analysis  of  concentration  vs  response 
data  based  on  a  logarithmic-logistic  function  for  both  mono-  and  dl-chotomous 
dose- response  curves. 

3.  Development  of  a  program  to  simulate  time  courses  of  drug  concentrations 
for  any  number  of  doses  and  route  of  administration. 

4.  Development  of  a  graphic  program  to  compute  area  under  the  curve  and  area 
under  the  first  moment  of  the  curve  for  drug  concentration  vs  time  data. 

5.  Development  of  a  graphic  program  to  compute  the  urinary  excretion  rates 
of  drugs. 

6.  Development  of  a  program  for  the  analysis  of  concentration  vs  effect  data 
using  probits. 

7.  Development  of  a  program  to  convert  pharmacokinetic  parameters  to  their 
exponential  forms. 

8.  Development  of  programs  for  the  analysis  and  simulation  of  non-linear 
pharmacokinetic  data. 

9.  Development  of  statistical  programs  useful  for  the  analysis  of  clinical 
pharmacological  data. 

10.  Work  on  a  program  to  help  design  drug  dosage  regimens. 

11.  Analysis  of  data  at  special  requests  of  COTR. 
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BACKGROUND 


The  process  of  drug  development  has  be«_n  both  complicated  and  facilitated 
.  y  the  trend  toward  early  application  of  the  methods  of  clinical 
pharmacokinetics.  Current  practice  and  existing  and  proposed  new  Food  and 
Drug  Administration  Regulations  demand  sophisticated  evaluation  of  drug 
L 1 oavl 1 abl 1 1 ty  and  descriptive  pharmacokinetics  In  Phase  I  clinical  studies. 
This  requires  development  of  assay  methods  for  new  drugs  and  their  metabolites 
and  methoos  for  evaluation  of  concentration  vs  time  data  to  obtain  relevant 
parameters  to  characterize  drug  behavior.  Of  similar  importance  is  knowledge 
of  characteristics  of  relationships  between  dose  or  concentrations  and 
pharmacologic  response.  The  Information  thus  obtained  from  pharamcokl neti c 
and  pharmacodynamic  studies  relates  directly  to  the  optimal  design  of  dosing 
schedules  of  new  drugs  1r,  man,  Including  individualization  of  therapy  due  to 
disease  processes  or  other  factors  which  may  affect  drug  behavior.  In  this 
regard  the  eventual  course  of  development  through  Phase  II  and  Phase  III 
clinical  studies  Is  rendered  less  empirical  and  ultimately  more  efficient  in 
both  time  and  cost,  by  minimizing  the  use  of  scarce  resources. 

Modern  computer  technology  has  greatly  enhanced  the  capability  of  these 
methods  and  made  possible  their  applications  to  clinical  pharmacokinetics.  A 
number  of  available  computer  programs  have  been  frequently  employed  for  this 
purpose.  Those  which  have  been  developed  for  use  specifically  In  clinical 
pharmacokinetics  are  based  on  com partmental  methods  of  analysis  and  yield 
estimates  of  parameters  associated  with  preselected  compartmental  models. 
While  useful  and  Informative,  they  lack  the  ease  of  use  and  cognitive  appeal 
of  direct  graphic  simulations  and  graphics-assisted  data  analysis.  It  was  In 
the  Interest  of  developing  general  purpose  programs  with  the  advantages  of 
computer  graphics  that  the  present  contract  was  Initially  pursued. 

Previous  experience  with  the  Tektronix  4051  and  4052  Microcomputer 
Graphics  System  Indicated  that  this  was  an  especially  suitable  microcomputer 
system  for  our  purpose.  While  similar  systems  are  now  available  from  a 
variety  of  sources.  It  Is  In  the  Interest  of  conformity  with  existing  systems 
In  the  U.S.  Army  Drug  Development  Program  at  the  Walter  Reed  Army  Institute  of 
Research  that  the  Tektronix  System  has  been  employed.  The  present  Tektronix 
4052  System  In  our  laboratory  and  for  which  the  programs  to  be  described  were 
developed  Is  Identical  to  that  now  In  use  at  Walter  Reed.  A  data 
communications  Interface  has  been  installed  to  facilitate  direct  transmission 
of  programs  and  data  between  these  two  facilities. 
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SUMMARY  CF  PROGRESS 


A.  Program  for  Linear  Pharmacokinetic  Data  Analysis  by  Non-Linear  Curve 

Fitting  Using  Compulsory  Graphics. 

The  foundation  upon  which  this  assembly  of  integrated  programs  was  built 
originated  several  years  ago  in  a  series  of  programs  developed  by  one  of  the 
investigators  (Dr.  Desjardins)  and  the  contract  monitor>  Dr.  Charles  Pamplin, 
in  the  Department  of  Pharmacology*  Division  of  Experimental  Therapeutics* 
Walter  Reed  Army  Institute  of  Research.  While  a  version  of  these  programs  had 
been  used  successfully  In  the  analysis  of  data  which  have  been  previously 
published*  they  had  not*  In  fact*  been  properly  verified  and  documented. 
These  programs  have  now  been  documented*  and  modified  and  expanded 
substantially  to  generalize  their  usefulness  to  the  greatest  extent  possible. 

In  the  course  of  this  work  a  concept*  not  entirely  original  but  new  in 
tne  extent  to  which  it  has  been  applied*  began  to  emerge  and  was  Incorporated 
In  the  program  package.  This  is  the  principle  of  separate  exponentials 
(Shand*  Desjardins*  Bjornsson,  Hammill  and  Pritchett:  The  method  of  separate 
exponentials:  A  simple  aid  to  devising  Intravenous  drug-loading  regimens. 
Clin.  Pharmacol.  Ther.*  29:542-547,  1981).  In  conjunction  with  the  unique 
compulsory  graphics  feature  of  the  program,  the  application  of  tills  principle 
has  resulted  in  the  assembly  of  programs  to  which  we  have  given  the  title* 
"Pharmacokinetic  Data  Analysis  by  the  Method  of  Separate  Exponentials  with 
Compulsory  Computer  Graphics". 

The  method  of  separate  exponentials  relies  on  Dost's  law  cf 
superimposition  for  resolving  the  exponents  Involved,  which  simply  states  that 
the  concentration  vs  time  curve  of  a  drug  is  made  up  of  the  sum  of  the 
exponential  terms  involved.  It  represents  a  departure  from  classical 
pharmacokinetics  as  it  has  evolved  over  the  past  two  decades  in  that  it  is  not 
directly  dependent  upon  compartmental  models.  It  provides*  Instead,  a 
descriptor  model  which  is  "statistically  sufficient"  to  define  the  behavior  of 
a  drug  with  respect  to  the  time  course  of  blood  (plasma)  concentrations.  The 
model  involves  a  nonlinear  regression  function  consisting  of  a  sum  of  the 
exponential  terms  which  des<  ri  be  the  total  time  course  of  the  drug 
concentrations  in  the  body.  The  model  provides  predicions  of  drug  behavior 
which  are  amenable  to  testing  by  suitable  clinical  pharmacokinetic  studies  and 
can  be  used  In  a  simulation  mode  to  optimize  drug  administration. 

The  equations  employed  for  each  route  of  administration  are  as  follows: 

1.  Intravenous  bolus  or  oral  administration 
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2.  Continuous  or  finite  infusion 
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3.  Oral  administration 
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i=2  (e  u ^  -e  ^1^)  force  through  0,  or 


n 

i  =  2  A^(e  ~  “  r^)  lag  time  x 

The  simplicity  of  these  equations  when  compared  with  compartment  model 
dependent  equations  with  Individual  1  ntercom partmental  rate  constants  Is 
Immediately  apparent.  They,  nevertheless,  when  statistically  valid,  provide 
estimates  of  relevant  parameters  and  serve  as  useful  predictor  models  of  drug 
behavior  in  a  linear  pharmacokinetic  system.  The  unique  feature  of  compulsory 
graphics  for  selecting  optimal  Initial  parameter  estimates  and  for  model 
selection  was  developed  In  the  original  programs  on  which  the  present  package 
Is  based  and  adapted  for  use  with  these  models.  Likewise,  a  series  of 
statistical  evaluations  for  model  selection  similar  to  that  originally  used 
has  been  adapted.  Modifications  include  the  use  of  a  combinatorial  function 
rather  than  a  table  of  numbers  for  the  sign  test  (Bury  KV:  Stastlcal  Models 
In  Applied  Science.  John  Wiley  and  Sons.  New  York  1975.  pp.  224-228)  and  the 
use  of  Akalke's  Information  crlterlon-AIC  (Yamaoka  K,  Nakagavla  T  and  Uno  T: 
Application  of  Akalke's  Information  Criterion  In  the  Evaluation  of  Linear 
Pharmacokinetic  Equations.  J.  Pharmacoki net.  Biopharm.,  6:165-175.  1978). 
The  latter  (AIC)  was  derived  from  the  maximum  likelihood  estimation  which  Is 
equivalent  to  estimating  the  parameters  so  as  to  minimize  "Kul  1  back-Lel  bl  er's 
mean  Information".  The  equation  with  the  minimum  AIC  Is  regarded  as  the  best 
representation  of  the  time  course  data.  Therefore,  this  method  Is  called 
"minimum  AIC  estimation"  (MAICE).  It  provides  a  statistically  valid  method 
for  estimation  of  the  number  of  exponential  terms  In  a  series  needed  to 
sufficiently  describe  the  fitted  data. 

The  parameter  estimation  routine  used  In  the  present  package  was  adapted 
from  Bailey  and  Homer  (Interactive  Parameter  Estimation,  Work  Unit  No.  M0098- 
PN. 01.0037,  Naval  Medical  Center.  Bethesda,  MD,  September  1976).  This 
program,  based  on  the  Marquardt  algorithm  for  least-squares  estimates  of  non¬ 
linear  parameters,  has  been  modified  slightly  to  provide  a  more  rapid 
convergence,  especially  with  the  optimal  Initial  parameter  estimates  provided 
by  "Compulsory  Graphics".  It  has  furthermore  been  verified  by  comparison  with 
the  BMDPAR  (Ralston,  1979,  UCLA)  program  for  pharmacokinetic  model  fitting. 
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B.  Program  for  Analysis  of  Concentration-Response  Kinetic  Data 

A  second  major  program,  first  developed  by  one  of  the  investigators  (Dr. 
Desjardins)  and  the  COTR  in  the  Department  of  Pharmacology,  Division  of 
Experimental  Therapeutics,  Walter  Reed  Army  Institute  of  Research,  involves 
describing  concentration-response  kinetics  by  nonlinear  curve  fitting.  In 
this  program  a  set  of  X,  Y  (concentrat 1  on- response)  data  are  fit  to  a 
logistic-logarithmic  function  which  was  parameterized  to  yield  a  direct 
estimate  of  the  IC-50  (50%  Inhibitory  concentration)  with  an  associated 
variance  estimate  from  which  an  approximate  95%  confidence  interval  Is 
calculated.  This  program,  which  had  been  used  in  the  analysis  of  aata  which 
were  subsequently  published,  had  not  been  properly  documented  and  verified. 

The  original  program  employed  the  graphics  capabilities  of  the  Tektronix 
4051  to  optimize  Initial  parameter  estimates  associated  with  a  hyperbolic 
tangent  form  of  the  logistic-logarithmic  function  as  follows: 

Y.  =  --y-  '  1  +  tanh  (  logC  -  loaX^)  }  +  L 


where 

Yj,X^=  Response  and  concentration  respectively  of  the  1th 
point  or  sample 
U  =  Response  at  X  =  0 
L  =  Response  at  X  =  ° 

C  =  IC-50 

r  =  A  slope  or  "scaling”  parameter 

This  Is  a  modification  of  a  similar  function  described  by  Emmons  In 
1942  and  Is  only  one  of  several  algebraically  equivalent  forms  of  the 
logistic-logarithmic  function,  all  of  which  describe  a  symmetrical  sigmoid 
relation  between  a  dependent  variable  Y  and  the  logarithm  of  an  independent 
variable  X. 

The  hyperbolic  tangent  form  of  the  equation  was  selected  originally 
because  It  could  be  written  to  Include  the  IC-50  as  one  of  the  parameters  of 
the  function.  Several  other  alternate  functions  were  subsequently  sought 
which  would  require  less  execution  time  and  would  also  be  flexible  with  regard 
to  Its  parameterization.  Three  general  types  of  equations  were  considered, 
all  of  which  provide  a  symmetrical  sigmoid  function  for  f(x,p),  1.e.»  the 
cumulative  normal  or  logistic  function,  the  arc  tangent  function,  and  the 
modified  logarithmic  logistic  function.  The  last  one  was  selected,  as 
fol lows : 


Yj ,  Xy  the  response  and  concentration  respectively  of  1^h 
samp! e. 


where 


Yc  -  the-  response  at  X  =  C 
V,  =  the  response  at  X  -  ■ 
f  =  (C  <  f<  1) 

ED^  =  the  concentration  associated  with  f  x  10C%  response 
S  =  a  scaling  (slope)  parameter 

This  form  of  the  logarithmic  logistic  equation  permits  a  direct 
estimation  of  any  response  point  with  its  associated  variance  estimate  from, 
which  a  confidence  interval  can  be  calculated.  It  also  offers  a  considerable 
degree  of  flexibility  and  ease  of  use  presently  not  available  In  other 
programs  for  dose-response  kinetic  data  analysis.  The  program  in  Its  present 
form,  following  input  of  X,Y  data  and  graphics-assisted  selection  of  Initial 
parameter  estimates  for  f  =  0.5,  will  optimize  these  parameters  by  nonlinear 
least  squares  curve  fitting.  Following  this,  the  data  will  automatically  be 
fitted  again  for  f  =  0.1,  f  =  0.9  and  f  =  0.95,  yielding  estimates  and 
variances  for  the  ED-10,  EC-50,  ED-90  and  ED-95.  With  minor  modification  of 
the  program  any  value  of  f  can  be  obtained  to  yield  an  estinate  for  ED-f.  The 
run  time  for  this  program  has  also  be  shortened  considerably  by  including 
exact  calculations  rather  than  numerical  approximations  of  the  partial 
derivatives. 

In  1980,  work  at  WRAIR  on  two  sets  of  data  pertaining  to  the  jLq  vitro 
concentration-response  kinetics  for  the  drug  PentostamR  against  promasti gotes 
of  Leishmania  donovani  showed  a  consistent  deviation  from  a  simple  sigmoid 
curve.  Work  on  this  phenomenon  resulted  In  the  development  of  the  option  In 
the  program  to  fit  the  data  to  a  "dichotomous  dose-response"  function. 
Conceptually,  the  model  assumes  that  some  proportion  of  the  organisms  (P) 
respond  to  the  drug  with  an  ID-50  =  Ct,  and  the  remaining  organisms  (1-P) 
respond  with  a  higher  ID-50  =  C2»  i.e.»  a  "dichotomous  dose-response". 
Mathematically,  the  model  thus  involves  two  dose-response  curves,  each  with 
its  ID-50,  slope  characteristics,  and  the  proportion  of  the  total  curve 
corresponding  to  that  part.  An  alternative  explanation  for  the  dichotomous 
dose- response  of  Leishmania  proma stl gotes  to  Pentostam  can  be  suggested. 
Namely,  It  Is  possible  that  the  drug  has  two  separate  mechanisms  of 
antll el  shmani al  activity  mediated  by  two  different  receptors  with  differing 
sensitivity  to  the  drug. 

A  new  program  was  therefore  developed,  based  on  the  existing  dose- 
response  program,  which  allows  the  analysis  of  data  exhibiting  a  dichotomous, 
or  two-step,  concentration-effect  relationship.  This  type  of  analysis  allows 
estimates  of  parameters  for  each  of  the  two  sections  of  the  dose-response 
curve,  1.e.»  two  ED-50's  are  estimated,  as  well  as  the  relative  fractions  of 
the  total  response  due  to  each  portion  of  the  curve.  This  program  has  proved 
valuable  In  analyzing  data  which  could  not  be  successfully  fitted  by  the 
conventional  dose-response  program. 


Plans  had  been  made  to  convert  the  dose-response  analysis  program  for  use 
on  a  TRS-80  computer.  This  was  not  possible,  however,  since  we  were  unable  to 
get  access  to  such  a  computer  for  the  required  amount  of  time  to  finish  the 
conversion,  and  it  was  not  felt  to  be  Justified  to  purchase  such  a  computer 
solely  for  this  project. 


C.  Program  to  Simulate  Time  Courses  of  Drug  Concentrations  after  Drug 
Administration,  for  Any  Number  of  Doses  and  Route  of  Administration 


A  comprehensive  and  versatile  pharmacokinetic  simulation  program  has  been 
developed  for  the  Tektronix  40o2  Microcomputer  Graphics  System  to  enable 
simulation  of  ccncentrati on  vs  time  curves  for  any  number  of  coses,  dosing 
intervals,  and  routes  of  admi ni stration,  and  any  changes  In  these  values.  The 
Tektronix  4050  Series  Microcomputer  Graphics  System  is  uniquely  well  suited 
for  such  a  program  because  of  its  exceptional  graphics  features.  The 
cevelopment  of  this  pharmacokinetic  simulation  program  was  a  logical  extension 
of  the  graphically-oriented  pharmacokinetic  data  fitting  program  previously 
developed.  Thus,  the  pharmacokinetic  simulation  program  can  serve  as  a  useful 
companion  tc  the  clinical  pharmacokinetics  program  package  developed  under 
this  contract  and  designed  for  use  in  the  U.S.  Army  Lrug  Cevelopment  Program. 

While  the  data  fitting  program  focuses  on  sophisticated  statistical 
analysis  cf  single  or  multiple  dose  pharmacokinetic  data  with  compulsory 
graphics,  inducing  simulation  of  the  best-fitted  concentration  vs  time  curve, 
it  has  been  designed  to  perform  pharmacokinetic  simulations  involving  selected 
kinetic  parameters.  Thus,  in  terms  of  use,  it  can  be  said  that  the  simulation 
4  program  takes  over  where  the  fitting  program  ends.  The  program  will  simulate 

predicted  drug  concentrations  for  any  given  dosing  regimen  when 
'  pharmacokinetic  data  are  available,  e.g.,  using  the  best-fitted 

pharmacokinetic  parameters  form  the  fitting  program,  or  pharmacokinetic 
parameters  from  the  literature. 

j  The  potential  uses  of  the  simulation  program  are  many  and  Include: 

1)  Design  of  drug  dosage  regimens.  This  Is  particularly  important  when 
designing  an  appropriate  chronic  dosage  regimen  of  drugs  or  drug  antagonist 
when  plasma  concentrations  have  to  be  maintained  In  a  certain  concentred  on 
range  or  above  some  minimal  effective  concentration,  either  for  prophylaxis  or 
treatment. 

2)  Illustration  of  the  effects  cf  intersubject  variability  In  drug 
deposition  on  drug  concentration  vs  time  curves  when  fixed  dosage  regimens 
are  utilized.  This  can  be  important  for  drugs  with  a  narrow  therapeutic 
index,  and  when  different  disease  states,  particularly  liver,  kidney,  and 

I  heart  disease,  may  affect  drug  disposition.  The  simulation  program  can  also 

,  be  used  to  predict  the  time  course  of  drug  concentrations  after  an  appropriate 

)  dosage  adjustments  has  been  made. 

[  3)  Optimization  of  blood  sample  collection  times.  This  Is  useful  for 

1  prospective  pharmacokinetics  studies  when  there  is  some  information  available 

I  on  the  pharmacokinetics  of  the  drug.  The  planned  sampling  times  can  then  be 

superimposed  on  the  expected  (based  on  best  available  information) 
concentration-time  curve. 

4)  Educational  uses.  This  includes  the  evaluation  of  the  effects  of 
changes  or  difference  in  individual  pharmacokinetic  parameters  on  the 
I  concenfati on  vs  time  profiles  for  different  dosing  regimens. 


The  user  of  this  program  has  several  options  available,  based  on  the  user 
cefined  keys.  Ifter  entering  the  pharmacokinetic  parameters  for  the  drug  of 


interest,  single  or  multiple  coses  may  be  entt  roc,  using  any  route,  amount, 
ane  time  of  cosing.  The  graph  nay  be  displayed  after  any  or  all  of  the  doses 
have  been  entered,  using  either  automatically  scaled  axes,  or  with  axis 
lengths  and  tic  intervals  set  by  the  user.  Logarithmic  scaling  of  the  Y-axis 
Is  also  available.  The  graphs  may  be  titled,  and  the  axes  may  be  labeled  at 
the  user's  discretion.  The  graph  may  be  drawn  on  the  screen  or  the  plotter, 
and  if  on  the  plotter,  the  program  will  pause  between  drawing  the  axes  and  the 
curve,  allowing  the  user  to  change  pen  colors  for  multiple  plots. 

The  output  of  this  program  consists  of  a  graph  of  the  predicted 
concentration  vs  time  curve,  and  a  printout  listing  the  pharmacokinetic 
parameters  used  for  the  simulation  and  the  doses  given,  including  amount, 
time,  and  route  of  administration. 

A  separate  simulation  program  has  been  written  for  the  simultaneous 
simulation  of  time  courses  of  drug  concentration  and  pharmacologic  effect. 
Conceptually,  this  program  involves  the  incorporation  of  the  log-logistic 
equation  into  the  linear  pharmacokinetics  simulation  program.  If 
concentrations  are  to  be  graphed,  they  may  be  shown  on  linear  cr  logarithmic 
scales,  while  the  effect  axis  is  always  linear  and  will  be  automatically 
scaled  by  the  program  (from  1%  to  100%  effect),  but  the  user  may  override  the 
choice  of  endpoints  and  tic  intervals  if  desired,  another  feature  which  has 
been  retained  from  the  original  simulation  program.  Another  useful  feature  of 
this  modified  program  is  the  pause  feature,  which  stops  the  pen  after  drawing 
the  axes  and  labels,  and  again  after  drawing  the  concentration  curve  and 
before  drawing  the  effect  curve,  allowing  the  changing  of  pen  colors  for 
different  sections  of  the  plot.  This  pharmacoki  r.eti  c/pbarmecody  nami  c 
simulation  program  should  prove  most  useful  in  designing  drug  dosage  regimens 
cr.  the  oasis  of  seme  desired  range  in  effectiveness. 

D.  Program  to  Compute  the  Area  under  the  Zero  Moment  (AUC)  and  the  First 

Moment  (AUMC)  of  the  Plasma  Concenratlon  vs  Time  Curve 

A  program  has  been  developed  for  the  Tektronix  4052  Microcomputer 
Graphics  System  to  compute  the  area  under  the  plasma  concentration  vs  time 
curve  (AUC),  and  also,  the  area  under  the  first  moment  of  the  plasma 
concentration  vs  time  curve  (AUMC).  These  two  parameters  are  cf  prime 
interest  in  pharmacokinetic  analysis  since  both  allow  non-com partmental  or 
model -1 ndependent  analysis.  The  parameter,  AUC,  can  be  used  to  calculate 
total  clearance  (Cl),  apparent  volume  of  distribution  by  area  (Vdarea)»  renal 
clearance  (Clp),  and  bloavall abil ity  (F),  using  conventional  expressions. 

The  area  under  the  first  moment  of  the  concentration  vs  time  curve  (AUMC) 
enables  the  calculation  cf  two  other  pharmacokinetic  parameters,  the  mean 
residence  time  (MRT)  and  the  apparent  volume  of  distribution  at  steady  state 
(Vdss).  These  parameters  are  calculated,  based  on  AUC,  AUMC,  dose  (D),  and 
duration  of  zero-order  drug  Input  (T),  e.g.,  by  continuous  i.v.  infusion,  as 
fol lows: 


voss  =  MRT  X  Cl 


which  In  the  case  of  i.v.  Injection  becomes: 

Vd  =  MRT  x  Cl  =  L_x_^'r!£ 
ss  ? 

Aur 


while  for  oral  administration*  MRT  has  to  be  computed  as: 


MRT  = 


AUMC  1_ 
AUC  K 


where  Ka  Is  the  first-order  absorption  rate  constant. 

The  MRT  represents  the  average  lifetime  expectancy  of  an  Individual 
molecule  within  the  body*  and  while  It  has  the  same  time  units  as  half-life* 
these  two  differ  significantly.  We  feel  that  this  parameter  will  be  of  great 
use  in  toxicology*  as  it  is  thought  to  express  the  exposure  of  the  body  or  a 
given  organ  to  a  given  compound.  The  other  parameter  derived  from  AUMC*  Vd5S, 
Is  a  volume  term  which  has  previously  only  been  calculated  after  compartmental 
analysis.  This  volume  is  usually  somewhat  smaller  than  Vdarea.  Thus,  the 
program  developea  to  computer  AUC  and  AUMC  allows  calcuations  of  Ci,  CIr* 
Vdarea’  and  Vdss'  addition  to  ^-1/2 

The  area  under  the  zero  moment  of  the  concentration  vs  time  curve  is 
calculated  using  the  trapezoidal  rule  over  the  range  of  observed  data  points, 
and  by  extrapolation  to  time  infinity  using  the  best-fitted  concentration  at 
the  last  time  point  and  the  terminal  elimination  rate  constant.  Values  for 
both  the  total,  trapezoidal  and  extrapolated  areas  are  listed. 

The  area  under  the  first  moment  of  the  concentration  vs  time  curve,  AUMC, 
Is  simlllarly  computed  by  summing  up  the  areas  under  the  first  moments  for 
each  trapezoid  from  time  0  to  time  t,  AUMC]*,  and  by  then  adding  the  area 
under  the  moment  curve  from  time  t  to  infinity*  AUMC]£  ,  which  Is  calculated 
from  the  equation: 

AUMC]  7  =  ~  *  ( t  +  l  ) 


where  C^  is  the  best-fitted  concentration  at  the  last  time  point*  t*  and 
the  terminal  rate  constant. 
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The  printout  generated  in  association  with  tie  analysis  lists  the  ca  t  a , 
parameters  of  the  final  phase  calculation  (intercept  ana  exponent)*  final 
half-life,  and  values  of  the  two  components  for  each  area,  plus  the  total  for 
both  areas. 

Through  the  user  definable  Keys  using  an  overlay,  the  user  has  the  option 
of  entering  data  directly  or  from  tape,  deleting,  changing  or  adding  data 
points  at  any  time,  giving  each  cata  set  a  title  which  will  be  printed  both  or, 
cop  of  the  graph  and  on  the  printout  (both  X-axis  and  Y-axis  can  also  be 
labeled),  and  deciding  which  data  points  to  use  for  computing  -.  Both  a  copy 
of  the  graph,  with  a  solid  line  indicating  which  data  points  were  used  to 
compute  1,  and  a  printout  will  be  generated. 


E.  Program  to  Compute  the  Urinary  Excretion  Rate  Constant  of  Drugs 

A  program  has  been  developed  for  the  Tektronix  4052  Microcomputer 
Graphics  System  to  compute  the  urinary  excretion  rate  of  drugs.  The  analysis 
yields  an  estimate  of  the  final  elimination  half-life  of  the  drug  based  on 
the  amount  of  drug  recovered  in  timed  urine  collections.  This  estimate  of 
rate  of  drug  el  1ml  nation  should  be  parallel  to  that  of  the  final  decline  In 
drug  concentration  In  plasma.  However,  unless  durations  of  urine  collections 
are  short  with  respect  to  the  elimlnaton  half-life,  and  the  different 
collection  periods  are  of  similar  lengths,  the  estimate  of  drug  elimination 
based  on  the  urinary  excretion  rate  method  tends  to  be  lower  than  that  based 
on  plasma  data.  The  reason  is  that  in  most  studies  durations  of  urine 
collections  Increase  with  time,  and  since  midpoint  values  for  urinary 
excretion  rate  are  calculated  for  each  collection  period,  the  elimination  will 
be  artifactually  decreased,  yielding  a  longer  final  half-life  value.  However, 
In  several  Instances,  urine  data  may  be  the  only  available  data. 

The  urinary  excretion  rate  is  computed  by  estimating  the  rates  of 
excretion  during  each  collection  period,  and  then  plotting  these  points  vs 
time  on  semi-log  axes.  The  values  on  the  Y-axis  are  therefore  /,  X  /  t  values 
for  each  of  the  midpoint  time  values.  The  elimination  rate  constant  is  then 
calculated  using  the  number  of  final  data  points  determined  by  the  user,  in  a 
manner  similar  to  that  used  In  the  area  program. 

The  printout  generated  In  association  with  the  analysis  lists  the  data, 
midpoints,  and  rate  of  excretion  for  each  collection  period,  and  the 
calculated  elimination  rate  constant  and  final  half-life. 

Through  the  user  definable  keys,  using  an  overlay,  the  user  has 
essentially  the  same  options  as  for  the  AUC/AUMC  program;  these  Include  direct 
data  entry  or  from  tape,  deleting,  changing,  or  adding  data,  entering  titles, 
and  generating  a  copy  of  the  graph  and  a  printout. 


F.  Progra*  for  Quantal  Dose-Response  (Probits)  Analysis 

The  probit  method  is  applicable  to  dose- response  analysis  when  only  a 
small  number  of  subjects  are  examined,  due  to  reasons  of  expense  or 
convenience  or  when  the  observed  response  cannot  be  precisely  measured.  This 
program  uses  a  probit  function  to  fit  the  observed  response  data  to  a 
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cumu !  cti  ve  normal  cistritution  function.  As  a  complement  to  our  other  ccsfc- 
respcr.se  program,  which  Is  to  sed  on  a  logistic  logarithmic  function,  this 
method  allows  the  user  to  analyze  data  from  exper ! nr.ents  in  which  the  frequency 
of  an  all-or-nothing  response,  such  as  death,  is  measured.  Interestingly, 
when  this  same  kind  of  frequency  data  Is  analyzed  by  the  logistic  cose- 
response  program,  the  results  obtained  are  quite  comparable  with  those  from 
the  ,.rGbit  program.  However,  ti  e  logistic  program  is  only  designed  to  be  used 
with  continuously  varying  effect  data,  while  the  probit  method  is  to  be  used 
for  discrete  responses.  Thus,  the  present  availability  of  programs  for  the 
analysis  of  both  graded  and  quanta!  responses  offers  means  of  analysing  all 
dose- response  data. 

The  method  used  for  analysing  quantal  dose- response  data  Involves 
converting  the  percentage  responding  to  each  dose  level  Into  probability 
units,  or  probits,  and  then  performing  a  linear  regression  of  protit  response 
vs  log  dose.  This  has  the  effect  of  fitting  a  cumulative  normal  probability 
curve  to  the  dose-response  data.  Due  to  the  nature  of  this  curve,  the  Initial 
fitting  is  usually  not  sufficiently  accurate  to  fully  describe  the  data.  To 
correct  this,  weights  are  assigned  to  each  data  point  based  on  its  position  on 
the  initial  regression  line,  and  a  new  weighted  regression  is  computed.  This 
process  Is  repeated  until  the  values  of  two  successive  interations  are  judged 
sufficiently  close.  Convergence  is  typically  reached  after  4-8  interations. 
The  resulting  line  corresponds  to  a  normal  probability  curve,  with  the  mean 
representing  the  ED-50  and  the  standard  deviation  giving  a  measure  of  the 
sensitivity  of  the  effect  to  changes  in  the  dose.  When  comparing  this  program 
to  the  logistic  pharmacodynamic  data  analysis  program,  we  obtain  very 
consistently  comparable  estimates  of  the  ED-50. 


G.  Program  for  Pharmacokinetic  Parameter  Conversion 

This  program  is  Intended  to  supplement  the  pharmacokinetic  simulation 
program.  A  difficulty  with  the  use  of  all  simulation  programs  is  the  widely 
varying  sets  of  parameters  which  are  available  in  the  published  literature. 
For  example,  it  is  not  uncommon  to  find  for  a  given  drug  that  the 
pharmacokinetic  parameters  provided  In  one  publication  are  co-efficients  and 
exponents,  in  the  second  the  model-dependent  compartment  parameters,  and  In 
the  third  a  selected  combination  of  both.  Since  our  simulation  program  Is 
designed  to  be  used  in  conjunction  with  the  sum-of-exponenti al  s  model  used  in 
the  linear  kinetics  analysis  program,  other  sets  of  parameters  must  be 
translated  into  this  format  in  order  to  be  useful  with  this  program  for 
pharmacokinetic  data  simulation. 

A  program  has  therefore  been  developed  to  automatically  convert  any 
sufficient  set  of  pharmacokinetic  parameters  Into  the  separate  exponential 
format  for  use  with  the  programs  already  developed  In  this  project.  The 
program  first  determines  whether  the  set  of  parameters  available  provides  a 
"complete  set"  and  If  so.  It  will  translate  them  Into  the  exponential  form. 
The  program  is  capable  of  converting  data  representing  one-  or  two-compartment 
pharmacokinetic  model  with  input  through  injection.  Infusion,  or  oral 
admi  ni  stration  into  a  single  standard  format  using  the  sum-of-exponenti al- 
terms  representation.  The  input  to  the  program  may  be  expressed  in  terms  of 
volumes  of  distribution,  half-lives,  time  of  peak  concentration,  total 
clearance,  rate  constants  of  absorption  or  elimination,  or  concentration  at 
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time  of  injection  or  maximal  absorption.  If  requested*  the  program  ear 
perform  unit  conversions  of  some  of  the  input  parameters  to  make  them 
consistent.  Not  only  will  this  translation  program  be  of  great  value  to  our 
simulation  program  but  also  to  all  Individuals  Interested  in  applied  clinical 
pharmacoki neti cs. 


H.  Programs  for  the  Analysis  and  Simulation  of  Non-linear  Pharmacokinetic 

Data 

The  programs  developed  earlier  for  curve  fitting  and  simulation  of 
pharmacokinetic  data  are  capable  of  fitting  and  simulating  time-concentration 
data  following  single  or  multiple  dose  administration  by  any  route  or  method 
of  administration  for  drugs  which  behave  according  to  linear-pharmacokinetic 
principles.  The  absorption*  distribution*  metabolism  and  elimination 
processes  are  all  assumed  to  be  first-order  and  their  net  effect  can  therefore 
be  described  by  a  sum  of  exponential  terms.  The  area  under  the  corresponding 
time  vs  concentration  curve  can  be  functionally  separated  according  to  the 
fractional  contribution  of  each  exponential  terms.  This*  Indeed*  Is  the 
foundation  of  the  concept  of  separate  exponentials. 

While  the  majority  of  drugs  in  current  use  and  new  developmental  drugs 
fit  the  concept  based  on  these  assumptions*  some  do  not.  Among  the  more 
common  examples  of  drugs  with  non-linear  pharmacokinetics  are  phenytoi  r,* 
sal Icyl 1c  aci d*  para-amlnobenzolc  acid,  ethyl  alcohol,  diphenhydramine  and 
sulfasalazine.  Nonlinearities  have  been  described  In  drug  absorption, 
distribution,  metabolism  and  elimination,  but  most  commonly  for  drug 
metabolism.  These  nonlinearities  imply  processes  which  are  not  first-order, 
but  rather  show  evidence  of  saturation.  The  net  effect  of  these  nonl 1 nearlties 
on  the  time  vs  concentration  curves  of  the  drugs  are  such  that  they  cannot  be 
accurately  described  by  a  simple  sum  of  exponential  terms.  They  require  a  more 
complex  model,  either  compartmental  or  physiological,  in  which  at  least  one  of 
the  functions  is  nonlinear.  This  function  Is,  In  turn,  most  commonly 
described  by  Ml chael  1 s-Menten  kinetics,  1  .e. : 


dC 

dt 


n 

i  =  1 
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The  program  developed  for  the  analysis  of  nonlinear  pharmacokinetic  data 
Is  a  modification  of  the  linear  pharmacokinetics  program  to  accommodate  sets 
of  differential  equations.  Instead  of  the  current  sum  of  exponentials.  A 
critical  component  of  it  has  therefore  to  be  a  rapid  and  convenient  method  for 
solving  the  differential  equations.  A  number  of  different  techniques  have  been 
evaluated,  including  Runge-Kutta  methods,  Hamming's  predictor-corrector 
system,  and  Bui  1 rsch-Stoer  rational  extrapolation.  None  of  these  has  proved 
to  be  Ideal,  In  the  sense  of  computing  solutions  rapidly,  and  being  easy  to 
use  and  compatible  with  the  other  graphical  and  statistical  components  of  the 
program.  The  existing  fully  functional  version  of  this  program  uses  a  fourth- 
order  Runge-Kutta  Integration  technique,  along  with  a  simple  means  of 
determining  step-size. 
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A  program  has  also  been  developed  for  simulating  time  courses  of  crucs 
behaving  In  a  nonlinear  pharmacokinetic  fashion  In  the  body.  This  program  Is 
a  modification  of  the  linear  pharmacokinetics  simulation  program,  but  Gue  to 
the  numerous  complexities  of  possible  nonlinear  steps  and  lack  of  available 
data,  the  present  program  has  only  been  written  for  the  one-compartment  model, 
with  Ml chael 1 s-Menten  type  elimination  and  first-order  absorption  when 
administered  orally.  Otherwise,  this  program  combines  the  flexibility  of  the 
existing  linear  pharmacokinetics  simulation  program  with  the  various  options 
available  In  the  linear  simulation  program,  such  as  mixed  modes  of  dosing  and 
user-sel ected  logarithmic  or  linear  scaling. 

In  the  present  version  of  the  programs,  the  user  does  not  have  to  enter 
the  differential  equations  describing  the  model,  but  can  now  choose  the 
appropriate  equation  from  a  menu  of  choices.  This  feature  takes  advantage  of 
the  ability  of  the  Tektronix  Basic  to  load  code  from  a  file  Into  an  executing 
program,  using  the  APPEND  statement.  The  nature  of  this  operation  makes  It 
very  easy  to  add  to  the  menu  of  models  by  adding  further  files  containing  the 
necessary  differential  equations,  and  then  editing  the  menu  to  reflect  the 
adcltlonal  choices  available.  Presently,  the  options  available  in  the 
standard  menu  of  models  Include  one-compartment  model  with  Michael  1 s-Menten 
elimination,  with  Input  from  IV  bolus,  IV  Infusion,  or  oral  first-order 
absorption.  These  models  are  the  most  commonly  used  for  nonlinear  kinetics, 
but  they  can  also  serve  as  a  template  for  the  development  of  further,  more 
complex  systems  of  differential  equations  for  these  programs. 

Together  with  the  already  completed  linear  pharmacokinetics  programs, 
these  programs  will  provide  a  comprehensive  package  of  programs  for  the 
analysis  and  slmmulatlon  of  pharmacokinetic  data  associated  with  virtually  any 
drug  of  Interest  to  the  U.S.  Army  Drug  Development  Program.  The  use  of 
compulsory  graphics  in  the  selection  of  model  selection  continue  to  remain  the 
unique  features  of  this  package  of  programs. 

I.  Statistical  Programs  for  Clinical  Pharmacological  Problems 

Work  on  the  statistical  program  package  began  with  an  emphasis  cn 
development  of  nonparametrlc  or  di  strl  btutlon-f  ree  techniques  for  comparing 
observations  made  under  two  or  more  courses  of  treatment.  These  tests  are 
applied  when  the  measurement  of  Interest  does  not  follow  a  normal  distribution 
or  when  the  number  of  observations  Is  so  small  that  the  nature  of  the 
statistical  distribution  cannot  be  determined. 

The  nonparametrlc  equivalents  of  the  two-sample  t-test  and  the  paired  t- 
test  are  the  Mann-Whltney  U-test  and  the  Wllcoxon  signed-rank  test, 
respectively.  Programs  for  performing  each  of  these  tests  have  been  prepared, 
with  both  being  designed  for  ease  of  use  with  respect  to  entering,  listing  and 
correcting  data.  In  each  of  these  programs,  the  user  Is  allowed  to  test 
either  a  one-tailed  or  a  two-tailed  hypothesis,  for  which  the  significance 
level  Is  then  calculated.  An  expression  for  the  exact  probabilities  has  been 
developed  and  Is  used,  thus  avoiding  the  need  for  tables  of  critical  values. 

The  nonparametrlc  methods  corresponding  to  1-way  and  2-way  analysis  of 
variance  are  the  Kruskal-Wal  1 1  s  1-way  ANOVA  and  Friedman's  2-way  ANOVA. 
Programs  for  these  tests  have  now  been  added  to  our  nonparametrlc  analysis  of 


variance  programs.  A  tic;  1 1  i  cna  1  1  y ,  a  program  has  been  written  for  the 
parametric  analysis  of  variance  with  multiple  comparisons  both  when  all 
possible  comparisons  between  groups  are  to  be  made  and  when  specific 
combinations  are  to  be  used.  The  multiple  comparison  methods  now  Involve  the 
methoc  developed  by  Schet'fe  and  the  Eonferrcni  method  of  adjusting 
significance  levels.  These  procedures  can  be  computed  directly  without 
requiring  the  use  of  tables  of  critical  values.  We  have  been  unable  to  find 
similar  direct  methods  to  use  with  the  methocs  of  Duncan  and  Dunnett.  Because 
of  this  lack,  these  procedures  are  still  not  available  as  an  automatic 
component  of  the  program,  although  the  user  can  easily  apply  either  of  these 
methocs  using  the  program  output  together  with  a  table  from  any  of  standard 
references . 


J.  Development  of  a  Program  Helpful  In  the  Design  of  Drug  Dosage  Regimens 

The  programs  developed  under  this  contract  for  the  analysis  and 
simulation  of  pharmacokinetic  and  pharmacodynamic  data,  1.e.»  drug 
concentration  vs  time  data  and  pharmacological  effect  vs  drug  concentration 
data,  respectively,  will  constitute  a  powerful  program  package  for  presenting 
clinical  and  basic  pharmacological  data  relevant  to  new  drug  development. 
These  programs,  however,  have  not  been  designed  to  address  directly  the 
Important  Issue  of  optimal  design  of  drug  dosage  regimens.  We  have  been 
Interested  In  developing  a  general  algorithm  for  designing  optimal  dosage 
regimens  for  drugs  for  which  both  pharmacokinetic  and  pharmacodynamic  data  are 
aval  1 abl e. 

Classically,  dosage  regimens  for  drugs  used  chronically  have  been 
determined  either  empirically  or  on  the  basis  of  pharmacokinetic  data  only.  A 
commonly  used  dogma  has  been  to  select  a  dosing  Interval,  ,  close  to  the 
elimination  half-life  of  the  drug,  t^/2»  1.e.» 


Method  I:  ~ 


This  approach  ensures  tnat  the  fluctuations  In  drug  concentrations 
between  doses  at  steady-state  will  neither  be  much  more  nor  much  less  than 
two-fold.  A  more  rational  method  would  be  to  select  on  the  basis  of  the 
desired  range  between  the  maximum  drug  concentration,  anc  the  minimum 

drug  concentration,  ,  over  a  dosing  Interval,  •  c  (  optimized  with 

respect  to  concentrations),  e.g.,  as  follows: 

In  (C  /C  ■  ) 

Method  II:  ■  ,  - — miji 

‘c  v 


where  K  Is  the  first-order  elimination  rate  constant.  Maintenance  doses,  D^, 
and  loading  doses,  D|_>  are  then  determined  based  on  pharmacokinetic  principles 
so  that  the  desired  drug  concentrations  aro  achieved.  Note  that  drug  dorage 


so  that  the  desired  true  concentrations  are  achieved.  Note  that  drug  cosage 
regimens  thus  designed  are  at  test  only  Indirectly  optimized  for  the  desired 
range  In  pharmacological  effect  between  doses. 


A  more  idealized  approach  would  be  to  base  the  selection  of  on  the 
basis  of  the  cesi red  range  in  pharmacological  effect  over  a  dosing  Interval, ' 

(  optimized  with  respect  to  pharmacologic  effect).  Such  relationships  have 
been  derived  by  one  of  us  (Dr.  Bjornsson).  In  the  simplest  case,  ‘  can  be 
determined  as  follows: 


(•E/m) 

Method  III:  -  =  - 


1.44  (.-E/m)  t1/2 


where  E  is  the  desired  range  In  the  pharmacological  effect  between  doses, 
e.g.»  between  60%  and  20%  effect,  m  is  the  slope  of  the  effect  vs  In  drug 
concentration  relationship,  and  K  and  t  w  2  are  as  defined  above.  Note  that 
this  approach  combines  information  obtained  from  both  pharmacokinetic  and 
pharmacodynamic  studies.  and  D|_  can  subsequently  be  determined  using 

pharmacokinetic  principles. 

Examination  of  Methods  II  and  III  reveals  that  for  any  given  value  of 
E,-e  Is  determined.  In  addition  to  K  or  two  by  the  slope  of  the  drug 
concentration  vs  response  curve.  Thus,  the  optimal  dosing  Interval  can  be 
considerably  longer  than  the  elimination  half-life  for  a  drug  with  a  shallow 
dose-response  curve,  and  vice  versa,  the  optimal  dosing  1 nterval  can  be 
considerably  shorter  than  the  elimination  half-life  for  a  drug  with  a  shallow 
dose-response  curve.  In  addition  to  these  pharmacokinetic  and  pharmacodynamic 
factors,  however,  other  factors  have  to  be  considered  in  the  design  of  the 
optimal  drug  dosage  regimen.  These  factors  Include  Interpatient 
pharmacokinetic  and  pharmacodynamic  variabilities,  therapeutic  Index,  toxic 
concentrations,  time-dependency  of  pharmacologic  effect,  and  convenience  or 
patient  acceptance  of  the  dosage  regimen. 

We  had  proposed  to  devise  an  algorithm  for  designing  optimal  drug  dosage 
regimens,  taking  Into  account  the  various  kinetic,  dynamic,  and  patient 
factors.  In  principle  such  an  algorithm  can  be  relatively  simple,  as  It 
Involves  simple  optimization  procedures  based  on  Methods  1 1 / 1 1 1 .  We  had 
anticipated  to  Identify  published  pharmacokinetic/pharmacodynamic  data  In  the 
public  domain  which  could  be  used  to  help  design  and  Illustrate  the  approach. 
However,  we  were  unable  to  Identify  a  useful  data  set,  and  this  project  was 
therefore  not  pursued  at  the  present  time.  This  Is  justified  as  more 
theoretical  work  will  now  need  to  be  done  on  this  approach  In  order  to 
overcome  the  lack  of  useful  data.  We  plan  to  continue  our  work  or  this 
exciting  project  In  the  near  future. 
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K.  Special  Requests  for  Data  Analysis 


In  addition  to  the  work  outlined  above  on  both  continuing  Improvements  of 
existing  programs  and  developments  of  new  programs,  several  special  requests 
for  data  analysis  were  made  by  the  contracting  Officer's  Technical 
Representative.  These  have  Included  the  following  Special  Reports,  which  have 
been  sent  to  the  COTR: 


1.  Data  Analysis  and  Report  on  a  Mefloquine  Blcaval  1  aMl  1  ty  Study  Comparing 
Formulations  E443  and  B514. 

2.  Data  Analysis  and  Report  on  a  Mefloquine  Bioaval  1 abll 1 ty  Study  Comparing 
Formulations  from  Lafayette  Pharmacal  and  Hoffman^-La  Roche. 

3.  Data  Analysis  and  Report  on  Blood  Chemistry  and  hematology  Studies  after 
Mefl oqul ne. 

4.  Data  Analysis  and  Report  on  a  Study  on  Methemoglobinemia  Caused  by 
Primaquine  and  Four  Other  Pharmacological  Agents. 


5. 


Data  Analysis  of  Ribavirin  Concentrations  In  Red  Blood  Cells  and  Plasma. 


